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CALCULATION OF LOW-PRESSURE INDICATOR DIAGRAMS. 

ByE. 0. Kbstblb. 

INTRODUCTION. 

It is well known that the recent remarkable advances in the design of high-speed internal- 
combustion engines have been made almost entirely by cut and try methods. Increased 
valve areas, increased compression ratios, better intake manifold design, better carburetion, 
etc., have enabled us to attain mean effective pressures that a few years ago were not dreamed 
of. Yet there is at present no accepted way of determining theoretically the relation between 
valve areas, stroke volume, speed, and volumetric efficiency. The question of the proper 
relative areas of inlet and exhaust valves is still moot. There is no way of fixing the proper 
valve timing for any given case, except the way of experiment, and there is little information 
available regarding the effect of improper valve settings on the volumetric efficiencies and 
pumping losses. In short, the designer dispenses almost entirely with theory. 

This condition of affairs is obviously due to the multiplicity of the factors which affect the 
perfoimance of internal-combustion engines. The large number of such factors not only makes 
theorizing extremely difficult, but it makes it exceedingly hard to interpret test results with 
any degree of certainty, so that we have little that is reliable in the way of data for checking 
up a theory once formulated. The increasing number of well-equipped experimental labora- 
tories, however, is sure to yield plenty of new and reliable experimental data, and just on this 
account it becomes important to develop a theory which shall guide the experimenter and at 
the same time aid in the interpretation of his results. 

It is hoped that the method of computing low pressure indicator cards outlined in this report 
will prove to be the first step in the development of such a theory. The writer particularly 
desires to emphasize the fact that this work is only a beginning, which can not at present be 
expected to yield results which are quantitatively reliable. 

The method of calculation here outlined is based on the following assumptions: 

(A) The mixture of air, burnt gas, and fuel in the cylinder may be treated as a perfect gas. 

(B) The kinetic energy developed by the incoming charge as it passes through the throat of 
the inlet valve is converted into heat by eddies as fast as it emerges from the valve, so that we 
may treat the gas in the cylinder as if it had a definite temperature and pressure. 

(C) The heat transmission between the gas and cylinder walls together with the heat 
absorption due to the vaporization of the fuel have a negligible effect on the pressure variations; 
that is, the processes under consideration may all be treated as adiabatic. 

(D) The rate of air flow through the intake and exhaust valves may be calculated from the 
manifold pressures, the cylinder pressure, and the valve clear openings by the ordinary formulas 
for steady flow through an orifice with the aid of determinable coefficients of effiux. 

These hypotheses are all to be regarded as rough approximations to the truth. The writer 
is of the opinion that the phenomena ignored through the introduction of these assumptions are 
of secondary importance in determining the pressure variations in the cylinder during the suction 
and exhaust strokes. In the last analysis, however, the justification of this opinion must come 
through a comparison of the results obtained from this theory with the results of experiment. 

Assumptions B and D are apparently indispensable to the present method of attack on the 
problem. A and C, on the other hand, constitute convenient simplifications suitable to the 
present preliminary analysis, but not necessarily essential to the method. In connection with 
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hypothesis 0 it should be observed that even though the charge in the cylinder absorb a good 
deal more heat from, the walls during the suction stroke than is required to complete the 
vaporization of the fuel, the increase in the heat content of the charge will not be accompanied 
by a correspondingly larger final pressure. This is because a slight increase in pressure due to a 
rise in temperature immediately reduces the rate of flow through the inlet valves. 

In order to apply the method it is necessary that pressures in the intake and exhaust 
manifolds at points adjacent to the valves and the effective clear valve openings shall be known 
functions of the time. By the effective clear opening of a valve is meant the product of the 
minimum sectional area of the passage through the valve and the coefficient of efflux. 

The use of the method is obviously restricted by the above requirements, since the manifold 
pressures are usually not known. In the case of a motor exhausting directly to atmosphere, 
however, the manifold pressure becomes the barometric pressure. Furthermore, certain types 
of multicylinder inlet manifolds have a nearly constant pressure, while a motor which draws 
air directly from the atmosphere and has its fuel injected directly into the cylinder will have 
atmospheric intake pressure, so that in these special cases the method is obviously applicable. 
In general the fluctuations in the manifold pressure may be expected to be small compared with 
those in the cylinder, and it is believed that the result of applying the method by replacing the 
actual varying pressure by its average value will be instructive, though quantitatively in error. 

The coefficient of efflux may be expected to depend on the shape of the valve and port, 
the direction of flow, the ratio of the valve velocity to the velocity of the air through the valve, 
and the product of velocity of approach and the rate of valve closure. The coefficients for 
steady flow may be obtained from such experiments as those described in Report No. 24, Air 
Flow Through Poppet Valves, by Lewis and Nutting. The coefficients applicable in cases of 
intermittent flow might be obtained from a study of experimental light spring diagrams with 
the aid of the theory developed below. The writer is inclined from a priori considerations to 
the opinion that the coefficients for steady flow should not differ very greatly from those appli- 
cable to the intermittent flow obtained in practice provided that the pressures used in making 
the calculations are those actually existing at points close to the valves. (The pressures in 
the intake manifolds and induction pipes are different at different points as well as at different 
times.) This is by no means certain, however (Of. "Pressure Drop Through Poppet Valves," 
by 0. E. Lucke, Trans. A. S. M. E., vol. 27, 1905), and it is desirable that further experiments 
be made to determine the actual values of the coefficients of efflux in cases of rapidly varying 
flow. It should be pointed out in this connection that the errors involved in using the steady 
flow coefficients may be expected to increase with the ratio of the valve velocity to the velocity 
of the air through the valve and with product of the velocity of approach and the rate of valve 
closure. 

NOTATION. 

Let V t v= instantaneous volume of cylinder in cubic feet and cubic inches, respec- 
tively. 

S = specific volume in cubic feet per pound. 
P, p = absolute pressure in cylinder in pounds per square foot and pounds per 
square inch, respectively. 
P r , p'= absolute pressure of air approaching inlet valve in pounds per square foot 

and pounds per square inch, respectively. 
P*, p" = absolute pressure in exhaust passage in pounds per square foot and pounds 
per square inch, respectively. 
0= absolute Fahrenheit temperature of cylinder contents. 
6' = absolute Fahrenheit temperature of intake charge. 
Ai, di = clear ope nin g of inlet valves in square feet and square inches, respectively. 
Gi= coefficient of efflux of inlet valves. 
Fuf^GiAi, Oia t = effective clear opening of inlet valves in square feet and square inches, 
respectively. 

A„fa a = clear opening of exhaust valves in square feet and square inches, respec- 
tively.? 
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0,= coefficient of efflux of exhaust valves. 
F e , f e ^0 9 A a , C e a» = effective clear opening of exhaust valves in square feet and square inches, 
respectively. 

Mt =rate flow through inlet valves in pounds per second. 
7l/e=rate flow through exhaust valves in pounds per second. 
<7=mass of air and gas in cylinder in pounds, 
f =time in seconds. 

p = crank angle from head and dead center. 
R = gas constant. 

Or = specific heat at constant pressure in mechanical unite (foot-pounds per 
pound). 

<7 T = specific heat at constant volume. 

Po, V», 0o = assumed initial values of P, V, and 8. 

GENERAL ANALYSIS. 

Case I. — Burnt gases escaping through exhaust valves. 

Treating the cylinder contents as a perfect gas (Assumption A), we write 

PV=GR8 (1) 

C=>C,+R. (2) 

The specific heat of a gas at constant volume is a measure of the increase of internal energy 
of unit mass for each degree rise in temperature. The product of the specific heat at constant 
volume (assumed constant), the mass, and the absolute temperature is therefore equal to the 
total internal energy. Consequently the rate at which energy accumulates in the cylinder is 

jj(CvG8). Substituting from (1) into this expression, we obtain 

i^^' pv y^'ii^- dv » 

The rate at which external work is performed on the piston is P jj- The rate at which 

the escaping exhaust gas carries away energy is [CJ +PS]M e = Cp9M 9 . The rate at which the 
incoming charge brings energy in is [Oft' +P'S']M l =- Cvd'Mi. 

Neglecting the heat transmission to the cylinder walls (Assumption C), we equate the sum 
of the rate of accumulation of energy and the net rate of dissipation of energy to zero 

if It (PV) +p °^ M < ~ <&' M i-°' W 



or 



or 



2T + f F 37 + ^M e ~O8'M^0. ( 6 ) 

In virtue of (2), this becomes 

G, 1 dV. O v 1 dP. C p if. fifl'jf, 

c^ r TJt + v^c v 'p'di + v^c T ir — ~py °> < 6 > 

1 d , p y 1 dV y M A , 0,8' Mj 
In case the inlet valve is closed, M t vanishes and 5f e becomes—^- Then 

dl Ioge P= " -raT log ' V+y dl lo & G - ®> 

Hence 

p 0 \v G 0 ) \p 0 ej ' 
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This is the pressure temperature relation for adiabafcic expansion which might have been an- 
ticipated from our assumptions without the aid of analysis. 

Having determined the relation between the temperature and pressure in the cylinder, we 

can compute the value of jjj (tog a P) for any given values of P, V, and F e . Transforming 
equation (8), we obtain 

But 

^--J&--2?.¥CPAP"). (11) 

Where the function % which denotes the theoretical rate of flow through a well rounded orifice 
of unit area, has the form 

when 

and the form 
when 

The values of the constants and ^ in the above equations are 
If -p- is less than its critical value (10) becomes 

a log. p =irL"3r +_ P~ * J (16 > 

+ mjrj(ffi} (i7) 

p" 

If p is greater than its critical value, (10) takes the form 

Equations are to be evaluated by a graphical integration. 

Case II. — Freaft charge entering cylinder through inlet valves. 
In case the inlet valve is open and the exhaust valve is closed, equation (7) becomes 

3 - f l0g.P=-7y Tf +(7-l)-p T - . (19) 

- - WCP W*>} (20) 
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pi 

The ratio p will generally be less than its critical value and we may use the following approxi- 
mate expression for * (Cf., Sanford A. Moss, American Machinist, Sept. 20 and 27, 1906) : 

*-f'-*.£V(£)*-(£F' " 



Let 
Then 

(20) becomes 



^ P -f-(,- MW /^f(i-.)(¥) ' « 

Case III. — Burnt gas escaping through inlet valve. 

In modern high speed internal combustion engines the pressure inside the oylinder will 
generally be greater than the pressure in the inlet passage when the inlet valve begins to open. 
Consequently we have to deal with outflow as well as inflow through the inlet passages. On 
the other hand the flow through the exhaust passages is always in the one direction. If the 
periods of valve opening overlap, the outward flow may occur simultaneously through both 
valves. This case can be dealt with by the addition to the right-hand member of (18) of a 
term to take care of the outward flow through the inlet valve. 

Thus 

i^--^^V ^V(57? -(£y 

GRAPHICAL INTEGRATION. 

The above equations (17), (18), and (23) are of such a form that it is useless to attempt an 
analytical integration. They are of the general first order type. 

%=ffr,y) (25) 

An accurate method for the graphical integration of this type of equation is given in Runge's 
"Graphical Methods." For our purpose, however, it will do to use a simpler means of getting 
an approximate solution. 

The equation. (25) defines not one curve, but a whole family of ourves, one of whioh passes 
through every point in the region of definition off (x, y). f (%, y) is the slope at the point x, y 
of the particular curve which passes through x, y. We define the combination of the point 
x, y with a short straight line through that point having the slope/ (x, y) as a lineal element of 
differential equation (25). / (x, y) being continuous, the directions of these lineal elements 
vary continuously from point to point. Hence by plotting a large number of lineal elements 
we can estimate the direction of the integral curve passing through any chosen point (see plot 
4). It happens that we are interested only in the particular integral curve passing through the 
point logePo? to. Consequently labor is saved if we proceed as follows: First plot the lineal 
element through log, P 0 , to. Choose a slightly larger value of t, say and by a prolongation 
of the above lineal element to the line s^fu get an approximate value of the ordinate correspond- 
ing to this abscissa. Plot the lineal element through this new point and through two other 
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points having slightly greater and slightly smaller ordinates, respectively. With these lineal 
elements as a guide, it will be possible to sketch in by eye with fair accuracy the integral curve 
from its initial point on the line x=to to its intersection with the line 3?=^, and to get an approxi- 
mate value of the ordinate of its intersection with another suitably chosen vertical line, say 
ar«=f,. By a simple repetition of this process we can extend the integral curve to any desired, 
value of f with a minimnTn of labor. 

When the relation between the pressure in the cylinder and the time has been determined 
the mass of fresh charge drawn into the cylinder during the suction stroke can easily be calcu- 
lated by means of a second simple mechanical integration. We denote this charge by the 
symbol 6%. It is equal to the gross oharge flowing into the cylinder through the inlet valves 
minus the mass which flows out through these valves during the initial period of pressure equali- 
zation. Let t' denote the time at which the inlet valve begins to open. Let t" denote the time 
when the flow through the inlet valve reverses. Let t" ' denote the time of inlet valve closure. 
Then 

Gi= -JJkwVPAP'W+f/tm&'APtft. (26) 

These integrals can be evaluated by means of a plan i meter. 

One of the equations (17) or (18) is applicable at all times when the flow through the 
valves is outward, even though part of the flow is through the inlet valve, although, if gas is 
escaping from the cylinder into the inlet manifold, we must replace in (18) by (F t + F a ). It is 
not to be expected that there will ever be any inward flow through the exhaust valves, or an 
inward flow through the inlet valves when the exhaust valves are open. Thus Equations 
(17), (18), and (23) cover all cases which will arise in practice. 

INTRODUCTION OP NUMERICAL VALUES. 

Before making a definite application of the above equations, we introduce some rough 
assumptions regarding the properties of the gas in the cylinder. Normal exhaust gas from an 
engine using gasoline as fuel may be assumed to contain about 13 .6 per cent by volume of water 
vapor, 12.4 per cent of C0 1} and 74 per cent of N t . The percentages by weight are then 8.5 per 
cent (5,0), 19.1 per cent (CO,), and 72.4 per cent (iV,). The specific heats at 200° F. and at 
constant pressure may be taken to be 0.448 B. T. U. per pound, 0.202 B. T. U. per pound, 
and 0.244 B. T. U. per pound, respectively. The specific heat of the mixture in mechanical 
units computed from the above rough assumptions is 198 foot-pounds per pound. The gas 
constant of the mixture is 53,7 and the specific heat at constant volume works out to be 

0 T = O p -B - 198 - 54= 144 

The ratio of the specific heats is 

t = 144 =t1 ' 375 ' 

We take the above values of R and y as the basis for our first rough calculation. At the 
high temperatures prevailing during the exhaust stroke a smaller value of y would probably 
be more accurate, while it is probable that a slightly larger value would be better for the suction 
stroke. 

Introducing these values, we find that the critical pressure ratio is 0.5334. Further; 

K^= 0.526; ^=2.094; Kf = 1.094. 

Equation (17) reduces to 

|i„ g . P= _i^[f + „. 626x63 . 7J1 ..vr^y-]. 

It is useful to replace log, P by log 19 p, V by v, and F, by/ e . Then 
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2.3026 |log w 2> = -±^[§+0.526 X 53.7 X W^fJ^ 
Iflfttf f 2 ^ +338.8/. J- 



(27) 

Another useful change is to introduce the crank angle as the independent variable instead 
of the time. We denote the crank angle measured from the head end dead center by <p. We 
take the unit of measurement for this angle to be 10 degrees. Let N denote the crankshaft- 
speed in revolutions per minute. Then 

Tt —6d-*°- 6N 

d . n 0.596nZv ,338.8 , / pY> ,1!S 'l 
5 log- p—— Lft+oSy/- VHft-) J' ' (28) 

Let J. —piston area in square inches. 
r =■ crank pin radius in inches. 
Z= connecting rod length in inches. 

t) 0 = compression volume in cubic inches. 
Then 

v=v 0 +.4{r(l-cos <p) + l-i/P-f* sin^} =v 0 + A (l+r)-A[r cos tp+^P-r* sin* rf; (29) 

dv^ Ar sin p f r cos y ) 

3;*° 5.73 | 1 + VP-r» sin 1 pj* (30) 

Equations (28), (29), and (30) enable us to determine the relation between pressure and 
crank angle from the time when the exhaust valve opens until the pressure ratio has been 
reduced to its critical value. From the latter point till the opening of the inlet valve we must 
use an equation of the form (18). 

Introducing the value 1.373 for y and ohanging the independent variable from t to <p we 
derive the following equation from (18). 

d , 0.596[~<Zv , nnK J* r -/p"\°- tua //V'Y' 1 " /»"V" 8 I. 

Making the same changes in (23) that we have made in equations (17) and (18), we derive 
the following equation applicable during the suction Btroke. 

EFFLUX COEFFICIENTS. 

It will be assumed that efflux coefficients derived from steady flow experiments, such as 
those of Lewis and Nutting (Keport No. 24), are applicable to intermittent flow. As a first 
approximation, the effects of valve design and seat angle, as well as the effect of reversing 
the direction of flow, will be neglected. In other words, it will be assumed that the efflux 
coefficient of any inlet or exhaust valve is a function only of the ratio of the clear opening of 
the valve to the net area of the valve passage (i. e. the gross area minus the valve stem area). 

Hot 1 shows the relation between the above-mentioned ratio and the coefficient of efflux 
as computed from the experiments of Lewis and Nutting on li-mch valves. Whereas the 
coefficients of efflux reported by Lewis and Nutting are based on the rough assumption that 
the clear opening of the valve is equal to the product of the valve lif t and the portTiircum- 
ference, the values used in the curve on Plot 1 have been corrected by the substitution of clear 
valve areas computed in a more accurate manner by means of trigonometrical formulas. 
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APPLICATION TO TYPICAL MOTOR. 

Plot 2 shows the effective clear opening (product of actual clear opening by efflux coeffi- 
cient) for the exhaust and inlet valves of a typical engine plotted against the crank angle. 

Plot 3 shows the cylinder volume v and its derivative -J^- as functions of <p. 

We assume that pressure at the moment when the exhaust valve begins to open (<p =» 130°] 
is 78.0 pounds per square inch (abs.) and that the temperature is 3,000° on the absolute Fah- 
renheit scale. 

Let the exhaust be direct to atmosphere. Then 

2> o =78.0 
0 O = 3OOO 4 
t> 0 = 149.7 

Equations (28) and (31) become 

d, Q.mrdv , 16,720. . IS „"1 . . _ _ , nnv 

^log 10 2> T\^ + ~ ! sr^f jm i' P> 27 - 55 (38) 

d « 0.596 Y dv , 98,250 , //14.7V 1 ' 8 /14.7Y«H . " t9 . ' 

Let us further assume that the impact pressure in the passage leading to the inlet valve 
has the constant value 14.7 pounds per square inoh, and that the impact temperature of the 
inflowing charge is 70° F. These assumptions set to one side the temperature and pressure 
changes which occur in the inlet manifold. Then 

2>' = 14.7 
fl'=530° 

Equation (32), for the pressure variation during the suction Btroke, becomes 
a, 0.596 f<2v 37,150 ^ // 14.7-p T, ™J 14.7 -p \) 

The graphioal integration of these equations for the determination of the relation between 
log 10 p and <p is shown on plots 4, 5, and 6. 

Plot 7 shows the low-spring pressure-volume curves for 1,000 revolutions per minute and 
2,500 revolutions per minute derived from plots 4, 5, and 6. 

It will be observed that according to these diagrams the exhaust valve opens a little early 
at 1,000 revolutions per minute and a little late at 2,500 revolutions per minute. The inlet 
valve closing time is better adapted to 2,500 revolutions per minute than to the lower speed. 

CONCLUSION. 

In conclusion it should be pointed out that we have here developed a method of com- 
puting theoretical low-spring diagrams which is capable of further development and refinement 
not justified by the experimental data now available. At present the absolute values deduced 
from the theory can not be relied on. It is to be hoped, however, that the results will prove 
to be of qualitative value and that future experiments will show us how to so modify the 
theory that it will give quantitatively correct values. 

O 



